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Abstract 

We calculate the high-temperature two-loop effective potential using a general 't Hooft 
background gauge. The dependence on the gauge-fixing parameter £ is investigated. The 
effective coupling constant at the critical temperature 53 (T c ) 2 is decreased considerably 
compared to the one-loop result, independent of £. 
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There are strong indications that the electroweak standard theory predicts a first or- 
der phase transition at the electroweak scale |l|]-|§. Generating the baryon asymmetry 
of the universe at the electroweak phase transition is an exciting possibility. A better 
understanding of this phase transition is required, however, in order to clarify whether 
this is indeed the case. The electroweak phase transition cannot be treated completely 
by perturbative techniques. Problems are caused by infrared singularities of the symmet- 
ric phase, requiring the summation of infinite sets of diagrams. Lattice simulations take 
care of this automatically, and therefore are indispensable tools for the investigation of 
the electroweak phase transition However, they are not well suited for the study 

of important physical quantities like the sphaleron transition rate, or the rate of critical 
bubble formation. These quantities have been studied in quasiclassical approximation to 
one-loop order. Two- loop calculations should help to control to what extend the corre- 
sponding results are reliable, though of course genuine nonperturbative contributions to 
the potential in the infrared have to be taken into account differently. 

Both sphalerons and critical bubbles are static field configurations, i.e. they do not 
depend on the imaginary time variable r. Therefore, it becomes useful to integrate out 
the non-static Matsubara frequencies first, to some order in the loop expansion. This 
perturbative expansion should be reliable because the nonstatic Matsubara frequencies 
become heavy at high temperature. The longitudinal component of the gauge field Aq 
develops a Debye-mass proportional to gT and may be integrated out as well p0| . 

The resulting three-dimensional effective theory is of course non-local. Usually higher 
derivative terms are neglected in the spirit of the high temperature expansion. One also 
neglects the Weinberg mixing and considers the action of the three-dimensional SU(2)- 
Higgs model: 



where we have introduced dimensionless coordinates and fields 



(1) 



x^ — , $ -> v$, A^vA . (2) 
gv 



The scale v is left open for the moment. 

The effective 3-dimensional gauge coupling is defined as 



9z{T? = 9 — . (3) 
v 

The gauge coupling g has been scaled out of the covariant derivative and the field strength 
tensor. The high temperature effective potential is 

F h t($ t *) = A(($t$)2_ u 2 $ t^ . (4) 

For compactness of notation we use a rescaled "A" and u Vq" . They are temperature 
dependent constants which correspond as A = Xt/q 2 and Vq = {vq{T)/v) 2 to the one-loop 



quantities used in reference | H[| , respectively as A = X-i/g 2 anc ^ ^ v o — —mi/v 2 to the 
two-loop quantities used in ||. 

We divide the fields into a background and fluctuations 

1 

$ -> $ + 53$, $ = ^= ipe, $ = —= (crl + m a T a )e (5) 

v z 

A^A? + g 3 A?, i = 0, A? = a? . (6) 
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(The T a are the Pauli matrices.) In order to describe critical bubbles responsible for the 
onset of the electroweak phase transition it is sufficient to work with only one nonvanishing 
background component ip of 3> in an arbitrary but constant direction e. We consider this 
type of background only. 

Integrating out the fluctuating fields in the loop expansion generates an effective action 
to be used to find the saddle point solutions corresponding to sphalerons and critical 
bubbles. To higher loop order this expansion will break down for small values of the 
Higgs field ip. Near the broken minimum, however, the loop expansion is expected to 
work quite well. Whether the saddle point actions can be estimated reliably depends on 
how important different regions in field space are for the corresponding solutions. This 
question deserves further studies. 

In praxi it is not possible to calculate the full effective action, but one has to expand 
it in some way, cut off the expansion and calculate the coefficient in powers of derivatives 
of the field. This expansion must of course break down at small values of the field ip, 
because the derivative operator d has the mass dimension 1, which must be compensated 
by powers of p^ 1 . Indeed calculating the contribution of a higher term in the derivative 
expansion to the effective action of a quasiclassical configuration one finds a divergent 
result, except for the potential and the dupdup term. However for the latter kinetic term 
one finds a Z-factor [jD] which in one-loop order is very strongly gauge dependent and 
therefore even this term, if considered separately, is rather unphysical. 

What seems to be needed is some less gauge dependent (nonlocal) combination of 
kinetic terms. As we will demonstrate, the effective potential is less gauge dependent. In 
a strict expansion in g% its extrema are completely gauge-independent. It will play an 
essential role in the case of inhomogeneous field configurations as well. 

In order to integrate out the fluctuating fields one has to fix the gauge. We choose as 
gauge-fixing condition the 't Hooft background gauge 

T a = D{A) i A a i + l -£(&T a ® - #t r «$) = Q. a f - ^&ir a . (7) 

The resulting effective action T[ip] depends in general on the gauge- fixing. Physical quan- 
tities should, on the other hand, be independent of it. This is due to the fact that they are 
described by extrema of T[ip\. Kobes et.al. ]l2| showed, that the gauge-fixing dependence 
of the effective action can be written as 

sm = ^sxiv] > (8) 

where ip represents all kinds of fields and <5X[</5] is a functional of the fields which can be 
calculated from the gauge-fixing condition and from the generators of the gauge transfor- 
mation. One reads of immediately that the value of the effective action is gauge-fixing 
independent for solutions of the equations of motion, as it should. 

An often raised objection against the 't Hooft background gauges is (see e.g. p3[ ) that 
the field (p used in the gauge-fixing should be of another type than the background field 
and should therefore not be varied in calculating the equation of motion 

V Sip ) q, =tp Sip 
Nevertheless equation (|8|) holds even for these class of gauges with T[ip] = T[ip,ip = ip] as 



is explicitly shown in reference [12]. Similar statements, although less general, have been 



verified long time ago [14, 15]. 
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As a consistency check, we shall explicitly demonstrate the gauge-fixing independence 
(i.e. the ^-independence) of the value of the effective potential at its extrema. In order to 
achieve this one has to be careful to work consistently to a given order of g 2 . 

Expanding in terms of the fluctuating fields we obtain the action 

Sbt + I d 3 x (\-T a T a + Cfp) -» / d 3 x\ (10) 



2£ 



1 
1 

9z 



1 fl r 

-diipdiip + V^-ip' 



+ 



dupdicr + A {^p 2 - vf\ atp 



V Iff 



+ {di^ a f + ^ Qfr* 2 + A (tp 2 - v 2 )^j ttV - c a *d 2 c a + \& 2 c a *c a + af^diif ) + 
53 ( ^e a6c a>5 (^ a i " + ^ abc a> b di7r c - e abc d i c a *c b a'i - ^e°* c (pir a c b *c c 

+A (a 2 + 7r a Tr a ) atp + la? (^ff7r a - ff3 i7 r a ) + ^ff« + ^ffc a *c a ) + 

4 



g 3 [ -e e a^-a^ + -0^71- + 7 A(ff + vr vr ) +-(ff +7T7T ) 



One reads off the propagators and vertices. The gauge boson (af ) propagator may be 
written as 



The Higgs (ff), Goldstone (7r a ) and ghost (c a ) propagators are as usual with the masses 

ml = A(3(^ 2 - v 2 ) , m 2 Gs = \(<p 2 - v 2 ) + -^ 2 , (12) 

m 2 w = \ip 2 , 171% = ^ip 2 . (13) 

Note that there are no IR-divergences due to massless Goldstone-bosons at the broken 
minimum, for £ 7^ 0. In the background field formalism the fluctuation appear only in the 
inner lines while the external lines consist of background fields [l(| . Only the lPI-graphs 
contribute to the effective action. 

In the following we study the effective potential. Up to now it has been calculated to 
two-loop order in Landau gauge from the four-dimensional theory [17, 18| and from the 
high-temperature theory ||. In a recent work M. Laine fll| presented a calculation using 
a general covariant gauge (T a = diAf). 

The aim of our work is to calculate the potential in two-loop order using the 't Hooft 
background gauge with an arbitrary gauge- fixing parameter £, and to investigate the in- 
dependence. 

At one-loop one obtains the effective potential 

y}s = A(^ 4 - \vW) " ^(™h + 3m 3 Gs + (6 + 3£ 3 / 2 )m^ - 6m 3 gh ) . (14) 
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Figure 1: The Graphs contributing to the two-loop effective potential. 



The two-loop potential receives contributions from the graphs shown in figure [|. We 
use the MS-renormalization scheme. The sunset graphs consist of integrals with three 
denominators; the figure-eight graphs have two denominators. The latter can be calculated 
easily while the momenta in the numerators of the sunsets cause some trouble. We removed 
them by a procedure which is similar to the one used in reference ||. The remaining 
integrals are of the following form (d = 3 — 2e) 

nt \ *e f d d k f d d p 1 1 1 n 

G(mi ' m2) = " J (2^7 J^rWTT^WT^ = T^ mim2 (15) 

rr, , 4e f d d k f d d p 1 1 1 

H(m 1 ,m 2 ,m 3 ) = [i / TTTZw / 77^7 ,o , o o I _ 2 7Z , , _ 2 ( 16 ) 



(2ir) d J (2ir) d k 2 + m\ p 2 + m 2 . (k + p) 2 + m| 
1 7 + M4,)U ' fl-ln r + "' 2+m3 |jn7 i 



32vr 2 V2e ' v 7 16vr 2 v 2 V A* 

The two-loop potential is a rather long expression and therefore not displayed here. 
All the ingredients are given in the appendix. While single graphs have ^-dependent 
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divergences the overall divergence is ^-independent 

- 6A V - 3A, 2 + §^ 2 . (18) 

From equation (|i~7|) one sees that the /i-dependent part of the potential is proportional to 
the divergence and therefore ^-independent. Hence it is possible to treat the /i-dependence 
for one value of £ and the ^-dependence for one value of The former is usually discussed 
by means of the renormalization group. This has been done for the Landau gauge elsewhere 
H and will not be repeated here. Instead we are setting \i to be the value of the field in 
the broken minimum later on. 

As mentioned above the value of the action must be independent of the gauge-fixing on 
its extrema. In our case the value of the potential at its extrema should be ^-independent. 
In order to expand it systematically in orders of g\ one has to expand the value of the 
field at the minimum first 

¥W = A + &<P%L + 2 3 V2n + 0(£) ■ (19) 
Plugging this into the effective potential one gets 

^min = V eS (ip min ) = V<$ + gjV^l + 9lV^ n + O(gl) . (20) 

We have to distinguish two cases 

(i) One tree- level extremum is at = and will stay there (</?min = ^min = 0)- 



From equation <Mt) on gets: 



Vi°) = , (21) 



mm 



v£L = -^(-*«3) 8/a ■ (22) 



= ^2^o(-3-8A + 41n(^2 v /-A^jj . (23) 

V^ n and V^ n are real for < 0, i.e. above the tree-level roll-over temperature, and 
complex below. This is due to the fact that y ffiit i = is a maximum of Vht(<p) for Vq > 0. 
(ii) At this temperatures one should expand around the broken minimum <£W n = v o 



J 1 ) 

t mm 



32Avr 



(l + 2 + 2 5 / 2 A 3/2 ) , (24) 



= , (25) 

V® = 

mm 



-^L(3 + 2 7 / 2 A 3 / 2 ), 3 , (26) 

y( 2 ) = 3 y 2 (_ 3 + (11 + 421n(2/3))A - 2 9 / 2 A 3 / 2 + 24A 2 - 2 n / 2 A 5 / 2 - 128A 3 
1 1 10247HA V 



8A (l - 4A + 8A 2 ) In (l + V2\\ + 64A 3 In (3V2X 



+ 

+ 2A(-17 - 16A + 64A*) In j . (27) 

The values of vffi, (i = 1, 2, 3) are independent of £ as they should. Note that we do not 
have any IR-divergences showing up in covariant non-background gauges fiT| , 



5 



2.5 l-r 




0.02 0.04 0.06 0.08 0.1 0.12 

A 

Figure 2: g% vs. A at the critical temperature 



Equations ( p5[ ) - (^7|) show that the effective expansion parameter at small A is j" 3 ^ • 
Consequently it becomes of order 1 if one approaches the critical temperature from below. 
Going to the limit A — ► small therefore does not help in improving the convergence of the 
loop expansion close to the critical temperature. The overall i in Eq. ( p7| ) arises entirely 

from inserting y£j n into the tree- and one-loop potential. 

Although the expansion around <p^ m = vq is manifestly gauge-fixing independent, it 
is not useful close to the critical temperature. As soon as tp^ m tends towards zero the 
true position of the minimum may no longer be considered as being obtained as a small 
perturbation around Vq(T). In the following, we therefore do no longer insist on the g% 
expansion of (p m \ a - Consequently, we do no longer work consistently to a given order in 
and some £ dependence must be expected. A small £ dependence would be an indication 
of a reasonable convergence of the approximation. 

Both the one loop and the two-loop potential predict a first order phase transition. 
They have two local minima, which are degenerated at the respective critical temperature. 
Up to now the value of v used to rescale the field in equation (g), is arbitrary. It is an 
appropriate choice to take v and fi to be the value of the scalar field at the broken minimum. 
The asymmetric minimum of the rescaled field ip is then at tp a = 1 

V*s(<p = l) = . (28) 

In addition, we have the condition 

V eS & = 0) = V eS (<p = l) . (29) 

at the critical temperature. g% and Vq, the two parameters of the high temperature effective 
action (|l]), can be calculated from equations (|28| ) and (p9|). From this it follows that the 
two coupling constants A and g\ are not independent at the critical temperature. Since g\ 



6 




determines the sphaleron rate this relation is important to determine cosmological bounds 



on the Higgs mass [20 



In Figure |2[ g\ is taken at the corresponding one and two-loop critical temperature, 
respectively, and is plotted versus A for £ = 0, 1, 2. The gauge-fixing dependence is weak, 
both at one and two-loop. The inclusion of the two-loop contributions changes the mag- 
nitude of 53 (T c ) 2 by about 40%. This would reduce the sphaleron rate by many orders 
of magnitude. Therefore bounds on the Higgs mass from the wash-out of the baryon 
asymmetry may be less reliably than thought so far. 

The large corrections to g%(T c ) 2 are not caused by large corrections to tp m i n (T) /T 
at fixed temperature, but by the shift of T c . This is demonstrated in figure H (The in- 
dependence is again not significant her.) Crosses denote the one (x) and two (+) loop 
critical temperatures. Going from one to two- loop one essentially moves along an almost 
universal curve. 

The determination of the critical temperature from the perturbative potential is of 
course very questionable because the latter is unreliable for small <p values. Still, it is 
remarkable that the two-loop potential leads to a lowering of T c and an increase of g% 
almost independently of the gauge- fixing parameter £. 

The one and two-loop potential at the corresponding critical temperatures is plotted 
for £ = 0,1,2 in figure ||. One first notices that the bulge between the symmetric and the 
asymmetric minimum grows from one to two loop, which is in agreement with previous 
calculations using Landau gauge [17, Accordingly, the phase transition is stronger first 
order, as predicted by lattice calculations (7], B|. While V^\ip) depends strongly on £, 

(2) 

the shape of V^ s (tp) is only slightly ^-dependent. This indicates that the loop-expansion 
might asymptotically converge towards an <!; -independent effective potential if one uses 
't Hooft background gauges. 
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Figure 4: The one and two-loop effective potential for A = 0.12 and different values of the 
gauge-fixing parameter £ 



v(T), the field value of the broken minimum used for rescaling (equations (Q, ^)), is 
not a physical observable. It turns out to be ^-dependent. We have calculated the Z-factor 
for the Higgs-kinetic term at one loop 



di J 3 if m Gs — m^y m Gs — £ 3 / 2 m 3 



Z H {ip) = 1 + ^ ; + 



4vr [ m Gs + m w ' m^y \m Gs - m hs ~ £ m w 

10-13V£ + 9£ 1 f dm 2 w \ 2 2 1 /<9m|\ 2 



+ 



16 (1 + v 7 ^) mw \ ^ / 16 m lh V ^ 



1 1 ( dmlX 1 1 fdmly 1 
+ 16 m 3 Gs \ dip ) + 48 y dip J J ^ 

It has a strong ^-dependence for small (/^-values and becomes even negative at some range 
as already pointed out in reference |ll]] . Near the broken minimum however it behaves well 
but is still ^-dependent. One can discuss^, if this ^-dependence cancels the one of v (T) 
calculating the renormalized field value \f Z{v(T))v (T). We found that there is indeed 
some 45% reduction if the one-loop Z-factor with the two-loop g 2 and the two-loop v(T) 
is used. Calculating the temperature dependent PF-mass would also require the Z^-factor. 

In conclusion, we have demonstrated the applicability of the class of 't Hooft back- 
ground gauges (eq. |7|) for studies of the electroweak phase transition. The main advantage 
of this class of gauge-fixings is the absence of IR-divergences in the broken phase, which 
are caused by massless Goldstone bosons in the class of covariant gauges [T a = diAf) 
which is prominent in literature. The latter class has been used by M. Laine in a recent 
publication [19| to calculate the two loop potential. He found that the loop expansion 
converges even in the broken phase only for small values of £. This is essentially due to a 



3 We thank C. Wetterich for raising this question. 
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^-dependent infrared divergence which does not show up in 't Hooft background gauges. 
In our opinion this shows the superiority of these gauges. The severe problems of per- 
turbation theory in the symmetric phase caused by nonperturbative condensates can of 
course not be cured either. 

Besides the improved IR-behavior in the broken phase there are some technical advan- 
tages due to the absence of a mixed Goldstone-gauge boson propagator. If one restricts 
oneself to the 't Hooft-Feynman gauge (£ = 1), the gauge boson propagator turns out to 
be quite simple as well (cf. eq. (11)). Note that the background gauges are also used in 
the computation of high energy cross sections j2l| . 

We showed explicitly that the value of the effective potential at its minima is inde- 
pendent of the gauge-fixing parameter order by order if it is expanded consistently in 

2 

93- 

Close to the critical temperature this expansion breaks down. Here we worked with 
the full two-loop potential and used a rescaling procedure which is especially suited to the 
treatment of quasiclassical solutions like critical bubbles |ll| and sphalerons pOfl . This 
procedure is not gauge-fixing independent but the ^-dependence becomes substantially 
weaker from one to two- loop order in the case of the effective potential (figure ||). 
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Appendix 



The propagators and vertices are read off equation ([n]). We used the form of the gauge 
boson propagator given in equation (|Tl|). It is straightforward to write down the two-loop 
graphs shown in figure |. While the "figure 8" graphs can be evaluated easily the momenta 
in the numerators of the "sunset" graphs have to be removed first. 

We did this in two steps: 
(i) after use of momentum conservation two loop-momenta k and p are left. The mixed 
scalar products kp are removed using the following identities: 

2kp _ k 2 +p 2 + m 2 



(k + p) 2 + m 2 (k + p) 2 + m 2 

d " f> kpF(p 2 ,k 2 ) = (32) 



(2ir) d 

/^CMW.*"> = i/^W,**) (33) 

(ii) The remaining momenta in the numerators are removed using the identities: 

k 2 m 2 



k 2 + m 2 k 2 + m 2 

d d k f d d p k m p 



1 - ( 34 ) 
(35) 



(2ir) d J (2ir) d k 2 + m 2 

Formula (|35| ) holds only due to the cancelation of IR-singularities. This reduction proce- 
dure has been performed with FORM. The integrals left are of one of the types given in 
equations (15, fLGj). 

The single graphs are given below where the combinatorical factors are chosen in a 
way that the contribution to the two-loop potential is given by: 

93 fi § ure 8 ' s ~ 2 ^ ^ sunsets ( 36 ) 
The dimension is d = 3 — 2e. 
WWW = 

H(m w ,rn gh ,m gh ) (-S^mly + -mly) + H(m w , m w , m w ) ( -9m^ d + — 



+G(m gh , m gh ) (-- + - £ - 3^ - + 3^ ) + G(m w , m gh ) (^ + 6£- + 6£d- 12 £ 
+ G(rovp,mw) f— 7 - 3 - + 3d 



4 d 



WGhGh 



H(m gh ,m gh ,m gh ) (^-^^mlv^j + H(m w , m gh , m gh ) ^+6£ - ^ rn^ 
+ G(m gh , m gh ) N-^ - - A + G(m w ,m gh ) (-3) 



11 



WGsGs 



H(m w ,m Gs ,m Gs ) (^+^m^ -3m Gs ^j + G(m w ,m Gs ) (^+- 
+ G(m Gs ,m gh ) (+7}t) + G(m Gs ,m Gs 



3 -2,2 



GhGhGs = 
+H(m Gs , m gh , m gh ) ( ~ ,~ 

WW = 



G(m gh ,m gh ) (-^fl + i^j + GK,m gh )(+3^ + 3^-6^ 

n< A f 9 3 1 9 i 3 a\ 



WGs = 



G(m w ,m Gs ) (-^ + -dj + G(m Gs ,m gh ) ^+-^ 



GsGs = 



+G(m Gs ,m Gs ) (V^ A ) 



WGsH = 

BYrn m m ^ 3 ^ ™H , 3 4s , 3 ™H , 3 „2 3 2 

(m w , m Gs , m H ) I -5 + 4 < + 4 <; + 4 " 2 ™ G 



3 1 4 



+ # (m Gs , m H , m gh ) ( +| m Gs m| - j -L m Gs - | 
+ G(m w> m Gs ) (+| + l-L- m Gs - \ -L- m^j 
+ G(m w , m H ) (+1 - \ m Gs + \ ±- m^j 
+ G(m Gs , m gh ) ( +7 C - 7 —2~ m Gs + 7 — 2~ m H J + G(m GB , m H ) f 
+ G(m H , m gh ) (+^ + ^ -gs " ^ m|) 

GsGsH = 

H(m Gs , m Gs , m H ) (+6 A 2 i^ 2 ) 



GhGhH 



F(m H ,m g h,m gh ) (~^£ 2 V? 2 ) 
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WWH 



3 2 1 4 3 2 1 2,^2 



H(m w ,m w ,mn) \+—ip — m H --(^ — g- m H + - ip d- - 

+ H(m w ,m H ,m gh ) \^ 2 ^f + §£y> 2 - - + 

V 8 8 16 16 

+ F(m H , m gh , m gh ) ( - 1 £ v? 2 — g- + | £ V 2 + Jr ^ 2 — j- m& 
+ G(m gh , m gh ) (+ A ^ 2 ^ " | ^ ^ m|) 

^ / 3 2 1 2 3 2 1 \ 

+ G{m w ,m w )\--<p -^- ma + - V 

( 3^ 2 1 3 2 1 \ 



+ G(m H ,m gh ) I -— f c^ 2 — g- + — <^ 2 — g- j 
^ 16 mfy 16 m^y 

HHH = 

H(m K , m H , m H ) (+6 A 2 9? 2 ) 

WH = 

G(m w ,m H ) (--r + q cM + G(m H ,m gh ) ( 



GsH = 

G(m Gs ,m u ) (+^ 

HH = 

+G(m H ,m H ) ( +- A 



GhGs 



G(m Gs ,m gh ) 
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